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ABSTRACT 


High spatial resolution spectrophotometric observations made in the wavelength 
region XX 0.6 - 2.0 fxm are used to study the Jovian and Saturnian limb darkening 
at 0° phase angle. Limb darkening coefficients (k) of the Minnaert function are 
derived for the cloud layers of both planets. A value of k = 1.0 is found for 
Jupiter over the entire disk while values of between 0.75 and 0.90 are found for dif- 
ferent latitudes for Saturn. These data are used to derive geometric albedoes (G) for 
the various belts, zones, spots and regions observed on Jupitet and Saturn. These 

It 

values of G and k are in turn used to show that an isotropic scattering model is invalid 
for Jupiter and that at least an asymmetric scattering function, such as the Euler 
function [o^ (1 + a cos0) ], is needed to fit the Jovian data. The Jovian scattering 
function is found to generally vary between 0 . 960 (1 + 0.85 cos 0) and 0 . 994 ( 1 + 0.42 
cos0) as a function of wavelength and the feature observed. The Saturn geometric 
albedoes and values of k indicate that Euler's function fails to adequately model the 

scattering properties of the Saturnian clouds. As a result it is suggested that 

r •" 

simple scattering theory may not apply to the Saturn clouds or that they 
are better represented by a cumulus cloud model. 
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IR SPECTROPHOTOMETRY OF JUPITER & SATURN 


This paper reports the results obtained from a series of spectrophotometric 
observations of Jupiter and Saturn in the wavelength region XX 0.6 to 2.0pm. The 
program was carried out in order to obtain information on the spectral albedoes and 
limb darkening properties of these Jovian planets. The empirical data provides a 

rs/ 

basis for determining the wavelength dependent single scattering albedo (c^) and 
scattering function for the bloud layers of each belt, zone or spot observed on 
Jupiter. Models for the clouds of Saturn are also suggested. 

I. INSTRUMENTATION AND OBSERVATIONS 

The observations were carried out with the Kitt Peak Observatory 50 inch 
(127 cm) reflector and a 10 channel PbS spectrophotometer which acquires data at 
10 wavelengths simultaneously (Binder and Jones, 1969). The wavelengths observed 
and the passbands are given in Table 1 . The observations of Jupiter were made with 
a 2 ? 1 76 entrance aperture so that the geometric resolution was on average 1/16 that 
of the equatorial diameter (39 M to 47") of the planet. Due to seeing effects and 
guiding errors during the 1 minute integration periods used, the effective resolution 
was reduced to 3" to 4", which was still sufficient to isolate the major bands, zones, 
and the Great Red Spot . 

Jupiter was observed on May 11, 1971, April 7, 1972 and on June 25, 1972. The 
phase angles of Jupiter and these dates were 2?4 10?8 and 0?1 respectively. 

Observations of X Serpentis, a GOV object, were made to provide calibration data. 

ft 

Saturn was observed with a 1 . 84 entrance diaphragm and hence the geometric 
resolution was about 1/11 the equatorial diameter (21") of the planet. As in the 
case of Jupiter, this resolution was degraded to 2" to 3" by seeing effects and guiding 
errors during the 2 minute integration periods used to make the observations . 

Saturn was observed on the 25, 27, and 28 of November, 1972 when the phase 
angle varied between 0?2 and 0?4. The stars, < Gemorium (G 8 lb) and rj Pegasus 
(G 211 - HI + F?) provided calibration data for the Saturn observations. 
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TABLE 1 . PASSBAND CENTERS AND PASSBAND WIDTHS OF THE 
10 CHANNEL PbS SPECTROPHOTOMETER 


Detector No. 


Passband Center, 
Mm 


Passband Width, 
Mm 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


0.594 
0.640 
0.725 
0.823 
0.964 
1.077 
1.240 
1.442 
1.554 
1. 992 


0.029 
0.021 
0.034 
0.059 
0.084 
• ’ 0.093 
0.094 
0.084 
0.071 
0.074 


A. 





H. NUMERICAL REDUCTION 


The numerical reduction of the data was carried out using the Kitt Peak 
National Observatory computer facilities . The data reduction procedures used 
for reducing the Jupiter and Saturn data are similar to those used for Mars as 
described by Binder and Jones (1972). The programs (i.e. Jupiter/Saturn vs 
Mars) differ only in that the geometric computations take into account the oblate- 
ness of Jupiter and Saturn, while Mars is treated as a sphere . 

Briefly, the program was used to convert the observed x and y coordinates 
of each spot for terrestrial atmospheric dispersion and then to compute a series 
of zenographic (and Saturnographic ) parameters [ latitude, longitude, cosine of 
the angle of incidence (i), cosine of the angle of emission (e), and reference air 
mass (AM)]for each point and to correct the photometric data for extinction, etc. 
In addition, the stellar data were corrected for wavelength dependent and seeing 
light losses since they were observed with 474 entrance diaphragm (the largest 
that can be used with the 10 channel system); the procedures followed are dis- 
cussed by Young (1970). Also, while the stars observed are solar or near solar 
types, corrections to their observed colors were necessary to account for the 
differences between them and the sun. The data of Johnson (1965) and Johnson 
et al. (1966) were used to make these corrections. 

m. APPLICABILITY OF THE MINNAERT FUNCTION FOR 
' LIMB DARKENING STUDIES 

In an earlier paper (Binder, 1972: hereafter referred to as paper 1) showed 
that the Minnaert function could be used to describe the Jovian limb darkening in 
the 1 . 5 pm region and that , when the effects of the strong 1 . 5 pm CH^ and NH^ 

absorptions were accounted for, the cloud layer appeared to behave like a 
Lambertian reflector. As a result of this work, we were interested in testing 
the general applicability of the Minnaert function to both Jupiter and Saturn over 
a wide range of wavelengths, to obtain a more precise determination of the limb 
darkening of the cloud layers, and to determine if the Minnaert function could 
also be used to compare theoretical scattering models with the observational data. 


The function has the form 


B = B cos^i cos^e * (1) 

so v ' 

where B is the brightness of the surface element, B the brightness of a surface 

8 k k-1 0 

element at cos i cos e = 1, i the angle of incidence, e the angle of emission, and 

k a dimensionless coefficient. Both B q and k are functions of the wavelength (X) 

and the phase angle (a). This equation can be written in the linear form, 

log (B . cos e) = log B + k log (cos i . cos e) (2) 

S 0 


so that k, which we shall refer to as the limb darkening coefficient, can be readily 
found from log-log plots of the data as demonstrated in paper 1 (see Figures 8 and 
9) and in Figures 1 and 2 of this paper using the new Jupiter and Saturn data. 

In addition to testing the general applicability of the Minneart function to fit 
empirical limb darkening data (see details in section IV below) we find that the 
function also satisfactorily represents the theoretical limb darkening for a variety 
of different Euler scattering models as computed by Horak and Little (1965), see 
Figure 3. As a result of these exercises, we find that the Minnaert function is an 
effective tool for representing and comparing empirical and theoretical limb 
darkening data. 

In addition, the Minnaert function is useful for deriving Harris’ limb darkening 
parameter x (1961), and the geometric albedo (G) of small areas on the disks of 
Jupiter and Saturn [both x and G are needed to determine single scattering albedoe 
(o^j) and scattering function of the clouds]. In the first case, Harris defines x as 



( 3 ) 


where B is the brightness of a surface element (B ) at the characteristic point, 
c s 

i.e. at cos i = 0.5, and B q is the brightness of the center of the disk . When a = Q, 


i.e. when cos i = cos e, 


„ 2k- 1 

B = B cos 1 
s o 


( 4 ) 


Thus, from (3) and (4) 

x = 2 (0.5) 2k_1 


(5) 


The geometric albedo (G) of a belt, zone or spot is derived from its brightness 
at the center of the disk, as computed from equation (1), as follows: 


G = /jBgds 
J 3 Fds 

where F is the incident solar flux per unit area. 

By equation (4) and (6), and assuming a circular disk, 

2tt 


( 6 ) 


G = 


r f r B cos 2k dir 
* / o o 

7 


2n Ff i r dr 


(7) 


where F is the flux per unit area . Now 
cos i = , 


/ 


( 8 ) 


by (7) and (8) 


B 0 fr d-r) 
F / r dr 


o k-1/2 . 
2 V ' dr 


Thus 


G = 


7 

±2 /Bq\ 

2 k + 1 \F ) 


(9) 


(10) 
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However B q /F is the normal albedo (N) of the point [ as derived from the high 
spatial resolution spectrophotometric data (Binder and Jones, 1972)] so the geo- 
metric albedo of any small point is 

G=2N (11) 

2 k+l 

For no limb darkening, i.e. k = 0.5, G = N, and for Lambertian limb darkening, 
i.e. k = 1.0, G = 2N/3 . 

Thus, via the Minnaert function and equations (5) and (11), spectrophotometric 
measurements of small areas of a planet can be readily converted into the quantities 
needed to study the scattering properties of the belts, zones and spots of the 
planets, even though the quantities themselves are not observed directly. 

IV. JOVIAN AND SATURNIAN LIMB DARKENING 

As indicated in Section III, and illustrated inpart in Figures 1 and 2, the limb 
darkening of all zones, belts, regions and spots of Jupiter and Saturn can be satis- 
factorily represented by the Minnaert function at all wavelengths.* Because of the 
axial symmetry of the markings of Jupiter and Saturn and because of the differences 
in center to limb and center to pole variations of the CH^ and (for Jupiter) NH^ 
absorptions (Hess, 1953; Munch and Younkin, 1964; Teifel 1969; Owen and Mason 
1969; Moroz and Cruikshank 1969; and Binder 1972) the E-W traces (i.e. traces 
taken along belts, zones and spots) must be treated separately from N-S traces 
taken along the central meridian. Table 2 and Figures 4a and 4b give the values of 
k for the various features observed and the averages of the N-S and E-W data res« 
pectively. 

As was found in the earlier study of the limb darkening of the 1 . 5 pm region 
(paper 1) and as apparent in Table 2, the limb darkening of the belts, zones, etc. 
are identical at a = 0 to within the uncertainty of the measurements, i.e. Jupiter 
has uniform E-W limb darkening over the entire disk. 

* Generally, however, the signal to noise for the 1.44 pm and 1.99 pm passbands is 
too poor for good values of k to be determined. This is the case since both of these 
passbands lie in very deep NH 3 and CH 4 absorption bands and telluric H 2 O bands . 


-TABLE 2 : LIMB DARKENING COEFFICIENTS FOR JUPITER 

AND SATURN 

Wavelength (jim) 0.59 0.64 0.72 0.82 0.96 1.08 1.24 
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The difference in the variations of k with wavelength for the N-S and E-W 
traces for Jupiter is readily explained as a result of the absorption variation data 
presented in paper 1. At almost all wavelengths, the limb darkening is due to the 
combined effects of the intrinsic limb darkening of the clouds and the center to limb 
absorption variations of the CH4 and NH 3 above the clouds. At shorter wavelengths, 
the contribution of the absorptions to the limb darkening is small (or zero); but 
with increasing wavelength, the absorptions become stronger so their contribution 
increases. The absorptions increase towards the limbs for N-S traces, thus the 
N-S limb darkening must also increase with respect to that of the clouds, i.e. 
k becomes larger, with increasing wavelengths. Similarly, the absorptions decrease 
towards the limbs for E-W traces, so the E-W limb darkening must decrease 
with respect to that of the clouds with increasing wavelengths. Thus, the observed 
values of k will converge, from opposite side, on the value of k for the clouds alone 
as shorter wavelengths are approached. Also from paper 1, we know that the absorp- 
tion variations in the IR are much stronger for the N-S traces then for the E-W 
traces, resulting in asymmetric, converging curves of k as seen in Figure 4a. 

This also implies that until the absorption becomes quite strong, i.e. for X > 1 pm, 
the limb darkening of the planet is essentially that of the clouds for E-W traces. 

Based on the above discussions and the fact that the N-S and E-W variations 

in k converge nicely to k = 1.00 at shorter wavelengths and the fact that measure- 

* 

ments of the limb darkening at or near the continuum level at 1.57 pm (paper 1), 
i.e. of the clouds alone, gives a value of k = 1.0 for N-S and E-W traces, we 
conclude that k = 1 or very nearly so (+ 0 . 01 ) for the clouds or haze layers over at 
least the wavelength range of XX 0 . 60 pm -1.60 pm . 

From Figure 4a and 4b it is clear that the limb darkening of Saturn is signifi- 
cantly different from that of Jupiter in several aspects. First, the N-S and E-W 
trace values of k do not converge at wavelengths less than 0.7 pm. Secondly, 

* The work reported in paper 1 was carried out at a spectral resolution of 300 A 
i.e. ~ 2.5 times the spectral resolution used in the 1.5 pm region for this work (see 
Table 1). As a result, one of the earlier passbands was nearly free from absorption, 
so the limb darkening of the clouds alone, at 1.5 pm, was determined from the earlier 
data. 


the values of k for wavelengths where CH^ absorptions are not important are less 
than 1 and there appears to be a latitude dependence of k for the E-W traces. Since 
Jupiter and Saturn seem so similar in most of their properties, it is quite significant 
that their limb darkening is found to be so different and thus the validity of the data 
has been checked to verify these findings . 

First, because of the comparatively smaller size of Saturn (20" vs 45" for 
Jupiter) and the presence of the bright rings, the effects of seeing and scattered 
light on the data are much more serious than for Jupiter. In order to determine the 
effects of seeing over the ball of the planet on the derived values of k, Dr. J. Pollack 
computed a series of seeing degraded limb darkening data for us . These computations 
show that for a worst case of 3" smear (due both to seeing and the finite diameter of 
our entrance aperature) our derived k's would be ~0.05 too high but within the error 
of our measurements. However, this correction to the derived k's would make the 
difference between the Jovian and Saturnian limb darkening even greater than observed. 
This same seeing smear would also bias the E-W limb darkening curves as a function 
of latitude, in that the higher latitude k's would be biased higher than the lower 
latitude k's; that is, in the direction opposite to the observed changes of k with 
latitude. Thus, the real differences in the k's for high and low latitudes would be 
greater than the observed differences if seeing smear of the ball were the only 
possible cause of error. 

The effects of the rings on the data are, qualitatively, as follows: Because of 
the geometry, the presence of the rings would tend to increase k for N-S traces, 
i.e. move it closer to Jovian limb darkening, and decrease k for E-W traces, more 
so for high latitudes than for low latitudes (this is the sense of the observed differences) . 
However, since at 1.24 pm the rings are brighter than the planet, the above effects 
would be maximized at this wavelength . Our data does not seem to show such an 
enhancement at 1.24 pm and so we conclude that the rings have also not seriously 
contaminated the data (i.e. such errors are on the order of or less than the observed 
uncertainties) . 
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Thus, the combined effects of the rings and seeing smear of the ball tend to 
increase k for N-S traces and tend to cancel each other on E-W traces. As a 
result, we conclude that the above noted differences between the Jovian and 
Saturnian limb darkening are real. 

i 

Based on Figure 4b and the above discussion, it is apparent that unlike the 
Jovian clouds, the Saturnian clouds do not behave like Lambertian reflecting layers. 
We conclude that the clouds have values of k in the range of 0.75 to 0.90 for the 
following reasons. First, at X < 0.7 pm, the derived k’s are essentially those for 
the clouds alone since there is little CH^ absorption at these wavelengths. Secondly, 
at X> 1.2 pm, the CH^ absorptions are completely or nearly saturated and so the 
brightness variations over the disk are those due mainly to the clouds. In fact, 
the 1.55 pm passband is completely clear of the CH^ absorption and so the k’s for 
X > 1.2 pm must also be essentially those for the clouds alone. At intermediate - 
wavelengths (0.7 pm < X < 1.2 pm) the CH^ absorption are not saturated and thus 
can effect the limb darkening, depending on the absorption variations over the disk 
(Teifel, 1969). Thus we can use the k's derived for X = 0. 54, 0.64, 1 .24 and 1 .55 pm, 
to derive the values of k for the clouds at low latitudes and high latitudes. In doing 
so, we note that for both sets of data, k at XX 1.24 - 1.55 pm is greater than k at 

4 

XX 0.59 - 0.65 pm. However, these differences are on the order of or less than the 
uncertainties of the observed k’s, so we do not consider them to be significant. 

_Thus, we assume k is approximately constant with wavelength and derive values of 
k = 0.85 + .05 for southern latitudes less than 35°and k = 0.80 + .05 for southern 
latitudes greater than 35 u for Saturn cloud limb darkening. 

V. SPECTRAL ALBEDOES OF JUPITER AND SATURN 

The photometric data obtained for each point on Jupiter and Saturn are converted 
into geometric albedoes via the observed values of k (Table 2) and the procedures 
discussed in Section HI. The geometric albedoes for the belts, zones, etc. observed 
are given in Figures 5 and 6 for Jupiter and Saturn respectively, as a function of 
latitude for all wavelengths observed except X = 1.44 pm and 1.99 pm for the 
reasons discussed in the note on page 6. In addition, Figures 7 and 8 give the 
geometric albedoes for representative features of Jupiter and Saturn as a function 
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of wavelength. These figures include, for comparison purposes, the geometric 
albedoes of the planets as derived from full disk photometry (Irvine et al. 1968 
and Newburn and Gulkis 1971). The photometric data are superimposed on 
schematic representation of the spectrum of each planet in order for the reader 
to assess the effects of the CH^ and NH^ absorption on the observed geometric 
albedoes. Discussion of the normalization of the spectrum to the photometric 
data is given in Section VI . 

There are four points to be made regarding the albedoes. First, there are 
significant differences between the two sets of Jupiter data represented in Figure 5, 
which are of course not surprising since the two sets are separated by 13 months 
of time. Previous long term observations of Jupiter (e.g. Irvine et al. 1968) show 
time dependent variations of 10% or greater in its geometric albedo. It is expected 
that individual features will show even greater time dependent variations than does 
the full disk. Thus, the data presented in Figures 5-8 has at least a 10% time 
dependent bias and that it is our opinion, that the current Jupiter data are on the 
lower side of the time average mean because of the general dusky appearance of 
the planet during our observing periods. 

Secondly, as is also known from color photography of Jupiter, the relative 

i 

contrast between the various features varies considerably with wavelength . For 
example in Figure 5 there are marked changes in the contrast between the South 
Tropical Zone (STrZ), the Equatorial Region [ the Equatorial Zone (EZ) and the 
North Equatorial Belt (NEB)] and the Great Red Spot (GRS) with increasing wave- 
length. These relative changes are in part due to the CH^ and NH^ variations 
over the features as discussed in paper 1 . 

Thirdly, the geometric albedoes given are not in general the albedoes of the 
continuum; the effects of the absorption due to CH^ and NH^ on the true albedoes 
are discussed in Section VI. 

* During both years, the EZ was nearly as dark as the NEB and so the entire 
equatorial region looked like a very broad, dark band which we refer to simply 
as the Equatorial Region (ER). 


Finally, because of the finite width of the entrance aperture, it was not always 
possible to measure a completely uniform area. Thus, the albedoes derived for 
small features such as the South Temperate Zone (STZ), South Temperate Belt 
(STB) and for the edges of the polar regions, which often contain faint zones and 
belts, are biased by their surroundings. 

/ 

In general the above argument also applied to the Saturn spectral albedoes, 
with one notable exception. As is shown in Figure 6, the contrast between the 
various features of Saturn does not appear to vary with wavelength. At all eight 
wavelengths Saturn has a relatively uniform albedo as a function of latitude except 
for a 10% dip at 30° south latitude . 

Also, the observed geometric albedoes of Saturn are higher than those of 
Jupiter. At wavelengths greater than 1 pm, this is in part due to the fact the NH^ 
absorptions are absent on Saturn and so the IR windows are more open for Saturn 
than for Jupiter (compare Figures 7 and 8). At shorter wavelengths, Saturn is 
10-20% brighter than Jupiter. This difference is probably, in part due to the 
time dependent variations mentioned above, but also may in part be time inde- 
pendent . 

We also note that the normalized spectral curves of Jupiter and Saturn 
(Figures 7 and 8) are very similar. Thus, we conclude that the color and composi- 
tion of the Jovian and Saturnian cloud layers are very nearly identical. 

VI. SCATTERING PROPERTIES OF THE JOVIAN CLOUD 
LAYERS 

The limb darkening data and geometric albedoes derived in the previous 
section can be used to derive new information regarding the nature and physical 
properties of the Jovian cloud layers. Following the discussions of Harris (1961) 
it is apparent that, while the Jovian cloud layer particles do not behave as isotropic 
scatters, an asymmetric scattering function, such as Euler's function, might be 
applicable for representing the scattering properties of Jupiter when limb darkening 
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data and geometric albedoes are both available. This scattering function has the 
form 

P (cos0) = at) (1 + a cos0) . (12) 

where o^, is the single scattering albedo, 0 is the scattering angle and a is a 

v 

dimensionless coefficient which varies between 1 and -1. When a is 1 the particle 
scatters predominately forward, when a = 0 the particle scatters isotropically, 
and when a = -1 the particle scatters predominately backwards. Harris has com- 
puted for a few cases the values of x and G, values which we need to define 05 and a 
(see Harris 1961, Table 23). Similarly, Horak and Little (1965) have computed 
the brightness over the disk of hypothetical planets for several additional cases, 
over a wide range of phase angles. However, these sets of theoretical data are 
insufficient for the purposes of this paper. Hence, additional values of x and G were 
computed for the range of 0 < op £ 1 and -1 ^ a ^ 1 for zero degrees phase angle 

i 

(a). Following Chandrasekhar (1950) the brightness at a point on a disk whose 
cloud particles scatter according to the Euler function is, at zero optical depth, 


I (M, <P, P q , <p q ) = 


p co F ( 
o o J 

4(P +/0 ( 


<Mm) ^(p o) - a 0 (p) e (p 0 ) 


a [(1-Ji 2 ) 1/2 H (1) (M)] [(1 t> 0 2 ) 1/2 H (1) op ] cos (n> 0 - ») } 


(13) 


where, p is cos i, p q is cos e, 0 the azmuthial angle of the point from the sub-solar 
point, < 0 o the azmuthial angle of the point from the sub-earth point; ^(p), 0 (p Q ), 0 (p), 
0 (p Q ) H^ (p) and H^(p q ) are tabulated in Chandrasekhar (1950). 

However, at a = 0, p = p Q and cos (<P Q - 0 ) = -1, bo equation 13 becomes 

! o 
0 2 (p) - a 02(p)- a (1-p 2 ) H (1) (tf 

. 



(14) 


Thus, by (3) and (14) x is 


x . = 21 (0-S) = 2 {* 2 (0.5)~ a e2 (0.5) ' °- 75 aH<1) (0,5)} ( 15) 

^o.i) * 2 d - » e 2 (1) 

Equation (5) can then be used to express the calculated values of x in terms of 
k, the Minnaert limb darkening coefficient. Rewriting equation (5) we have 

k = 0.693 -ln(f) (16) 

1.386 

Also by equations (11) and (14) the geometric albedo is given by 



Using equations (15), (16) and 17), values of k and G were computed and are 
presented in Figures 9, from which values of a*) and a for each belt, zone or spot 
observed can be derived given the geometric albedoes of the clouds alone. 

As discussed in Section V, the albedoes given in Figures 6 and 7 for various 
areas on Jupiter are the observed geometric albedoes. These values are corrected 
for the NHg and CH^ absorption in the passband of the 10 channel spectrometer as 
follows. First, using spectrophotometric data on the IR spectrum at the center of 
the disk of Jupiter (Cruikshank and Binder 1968, 1969 and Woodman, private 
communication), a filter factor was derived which relates the observed intensity at 
each each wavelength (except 1 . 44 and 1 . 99 pm) to the maximum intensity of the 
window in which the passband lies . For example (see Figure 7) the maximum of 
the peak of the 1.5 pm window of Jupiter is found to be brighter than the observed 
intensity as defined by the 1.55 pm passband by 40%, or the maximum at 0.82 pm is 
19% brighter than that observed by the 0.825 pm passband, when the effects of the 


M. 


I 


absorptions are considered. Using these filter factors (see Table 3) and the 
photometric data fiven in Figure 5 were used to define the normalized Jovian 
spectrophotometric curve as shown in Figure 7. 

Secondly, the peak of the windows of Jupiter beyond 1.1 pm do not lie at the 
continuum level and so, using laboratory data (Cruikshank and Binder 1968) an 
additional coefficient was derived (see Table 3) which, when multiplied by the 
observed geometric albedo and the filter factor, gives the geometric albedo of 
the cloud layer. It is recognized that these filter factors, estimates of continuum 
level factors are difficult to make, especially at longer wavelength. where the 
absorptions are very strong. Thus, the final values derived for the geometric 
albedoes of the clouds are most probably in error by a few percent at short 
wavelengths and up to 10% or perhaps more at longer wavelengths . 

With these cautionary remarks in mind, values of G for the cloud layer of 
the various belts, zones or the Red Spot can be derived using the data in Table 3 
and the values of the observed G in Figure 5. Using these final values of G and 
a value of k = 1 for Jupiter, and a can be determined for the Euler function for 
each feature from Figure 9. 

Figures 10a and 10b give the values derived by this procedure for the STrZ 
and the dark EZ of Jupiter. From Figures 10a and 10b, it is apparent that the 
Jovian cloud layers show a strong wavelength dependency in their scattering 
properties. First, after an initial rise between 0.60 and 0.65 pm, drops off 

rsj 

rapidly with increasing wavelength to 1.2 pm where then apparently begins to 
increase moderately again. This behavior is not that expected for a uniform 
population of white scatters, since in this case would show a continual increase 
with increasing wavelength. However, since the cloud particles are most probably 
NHg crystals, it is likely that the decrease of pf, with wavelength beyond 0.65 pm 
is due to the solid state absorptions of the NH^ crystals, absorption which should 
increase in strength in the regions beyond 0.65 pm and hence darken the particles. 
Alternatively, it is possible that at the shorter wavelength, we are measuring the 
scattering properties of an upper cloud layer (Clements and Tomasko, private 
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TABLE 3 


FILTER AND CONTINUUM LEVEL FACTORS 



Jupiter 

Jupiter 

Saturn 

AOzm) 

Filter Factor 

Continuum Level 
Factor 

Filter Factor 

0.59 

1.00 

1.00 

1.00 

0.64 

1.03 

1.00 

1.06 

0.72 

1.29 

1.00 

1.54 

0.82 

1.19 

1.00 

1.29 

0.96 

1.45 

1.00 

1.55 

1.08 

1.30 

1.00 

1.15 

1.24 

1.51 

1.13 

1.31 

1.55 

1.40 

1.22 

1.00 
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communication). With increasing wavelength, the scattering efficiency of this 
upper layer may become increasingly smaller and so the light reaches a second 
cloud layer, whose is less than that of the upper layer. Thus we might be 
seeing the combined effects of two different populations of scatters at different 
levels in the atmosphere . 

Figure 10b shows that, at the shorter wavelength, the cloud layers are 
moderately forward scattering and that as longer wavelengths are reached, the 
scattering becomes more and more forward--. Using the data presented in Figures 
10a and 10b, we find that 

0.960 (1 + 0.85 cos 0) £ p (cos 0) £ 0.994 (1 + 0.42 cos 0) (18) 

for the general range of G over all observed wavelengths . The maximum and 
minimum values, based on our present data are 

p (cos 0) = 0.997 (1 + 0.30 cos 0) (19) 

r 

for X =0.64 pm in the STZ and 

p (cos 0) = 0.945 (1 + 0.94 cos 0) (20) 

for X = 1.55 pm in the NPR. If it is assumed that the maximum time dependent geo- 
metric albedoes are 10% higher than our value, then at X = 0.64 pm,one would have 
Gmax ~ 0 . 64 and from Figure 9 

p (cos 0) = 0.998 (1 + 0.20 cos 0) (21) 

Thus, even for the most favorable case, we find that a conservative, isotropic 
scattering model is not valid far Jupiter .and hence at least an asymmetric scattering 
function, like Euler’s function, is needed. While the more general form 

p (cos 0) = 1 + Wj P 2 (cos 6) + P 2 (cos 0) (22) 

may ultimately be required to fit Jovian scattering data (i.e. photometric 
data taken from spacecrafts at large phase angles or absorption line profiles taken 
at the extreme edge of the disk), the two parameter Euler function adequately fits 
the available data. 
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VH. SCATTERING PROPERTIES OF THE SATURNIAN CLOUD 
LAYERS 

In general, the Saturn data were treated in the same way as were the Jupiter- 
data in that the observed geometric albedoes were converted into continuum albedoes 
using filter factors (as given in Table 3) derived from IR spectral data of Saturn 
(Johnson, 1970 and Younkin, private communication). These data were used to 
derive the normalized Saturn spectral reflectivity curve given in Figure 8. Figure 
10 gives the ranges of and a for the Euler function, based on derived ranges of 
k and G for Saturn. 

From Figure 11 it is apparent that for k - 0.80 and 0.85 and 48% ^ G < 68%, 

(as derived from Figure 6 and Table 3) the Euler scattering function does not 
satisfactorily represent the scattering properties of Saturn. First, the function 
fails to give albedoes as high as those observed around 0.7 pm if, as is observed 
for high latitudes, k is less than about 0.85. This problem may not be fatal, 
because of the possible biases in both G and k. If Gmax were only ~ 60% and kmax 
~ 0.8, then the Euler function would satisfy all the daJta. However, at all observed 
wavelengths, -1 < a -0.3, so the function is predominately back scattering. 

This seems to be inconsistant with the general behavior of scattering particles. As 
a result, it is possible that a different function is needed to properly model the . 
Saturnian cloud layer scattering, or that simple scattering from uniform cloud 
layers does not apply in this case. For example, Schoenberg (1925) has shown 
that at a = 0 a Lambertian sphere whose surface is covered with small cones or 
hemispheres (which also reflect according to Lamberts Law) has values of 
0.7 <\ k ^ 0.85. The exact value of k depends on the shape of the proturbences and 
their relative spacings. Based on (1) the similarity of the derived values of k for 
Saturn's clouds and those derived by Schoenberg, and (2) the similarity of the 
geometry of cumulus clouds and to that of the cones and hemispheres used in the 
Schoenberg calculations, we suggest that Squires' cumulus cloud model for Jupiter(1957) 
may actually apply to Saturn rather than Jupiter. If this is the case, then the cloud 
particles of Saturn may individually behave like those of Jupiter [as is indicated in 
part by their similar reflectivity profiles (Section V)] but the macroscopic 
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"topography" of the Saturnian clouds cause the limb darkening of the two planets 
to differ. Also, the latitude dependence of k for Saturn may be due to changes in 
the geometry and spacing of the cumulus clouds with increasing latitude. 

VIII. SUMMARY 

Spectrophotometric data taken over the wavelength range XX 0.6 - 2.0 /im 
are used to show the general applicability of the Minnaert function for representing 
Jovian and Saturnian limb darkening. These data also indicate that the observed 
wavelength dependent variations of the limb darkening of both planets is really 
explained as the combined effects of the intrinsic limb darkening of the clouds 
and the center to limb absorption variations of NH^ and/or CH^. Also the limb 
darkening of the clouds of Jupiter is significantly different from that of the 
Saturnian clouds. The value of k for the Jovian clouds is found to be 1.0 for all 
parts of disk, while k for the Saturnian clouds is found to be latitude dependent and 
to be in the range of 0.75 to 0.90. 

Geometric albedoes are derived for various features of Jupiter. These albedoes 
coupled with the limb darkening data are used to derive and a for the Euler 
scattering function. As a result, we show that isotropic scattering is not valid . 
for Jupiter, but that at least an asymmetric scattering function, such as Euler's, 
is needed to fit the valuable data. The derived values show that 05 decreases with 
increasing wavelength. This decrease may be due to the color of the NH^ crystals 
which probably make up the clouds . 

It is also shown that the Euler function does not adequately fit the Saturn data. ' 
As a result, it is suggested that Squires' (1957) cumulus cloud model may apply to 
Saturn. 

This research was sponsored by the National Aeronautics and 
Space Administration under contract NASW - 2363. 
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CAPTIONS 


Fig. 1. Plots of log (B . cos e) versus log (cos i . cos e) for representative E-W 

s 

and N-S traces over the disk of Jupiter. 

Fig. 2. Plots of log (B . cos e) versus log (cos i . cos e) for representative E-W 

s 

and N-S traces over the disk of Saturn. 

Fig. 3. Plots of log (B . cos e) versus log (cos i . cos e) for theoretical planets 

s 

with thick cloud layers whose particles scatter according to Euler’s function. The 
bottom three curves are for a = 0. The upper most curvte is for a = 10°. 

Fig. 4a - 4b. Plots of the observed values of k for Jupiter and Saturn respectively 
versus wavelength and wavenumber. In both figures, the filled circles represent 
the k's for N-S traces and the open symbols represent the k’s for E-W traces. 
Further, in Fig. 4b, the open circles are the k's derived for E-W traces, for 
southern latitudes less than 35° and the open squares on the k's derived for E-W 
traces for southern latitudes greater than 35°. 

Fig. 5. Plots of the geometric albedo versus zenographic latitude for the areas 
observed on the central meridian of Jupiter . The open symbols are for visually 
bright features and filled symbols are for visually dark features. The circles 

4 

(open and filled) are data taken on May 11, 1971 and the triangles (open and filled) 
are data taken on June 25, 1972. The features observed are identified at the top. 
The wavelength is given on each graph. 

Fig. 6. Plots of the geometric albedo versus Saturnographic latitude for the 
areas observed on the central meridian of Saturn. The open circles are data 
taken on November 25, 1972 and the filled circles are data taken on November 27, 
1972. The wavelength is given on each graph. 


Fig. 7. Plot of a geometric albedo versus wavelength for the STrZ of Jupiter. 

The curve is the normalized spectral reflectivity curve of Jupiter. The crosses 
are the new data presented in the paper (the width of the crossbar is equal to the 
passband width of that detector) . The filled circles are data presented by Newburn, Jr 
and Gulkis (1971) and the open squares of data presented by Irvine et al. (1968). 

Fig. 8 . Plot of geometric albedo versus wavelength for Saturn (average of all points). 
The curve is the normalized spectral reflectivity curve of Saturn. The crosses are 
the new data presented in this paper (the width of the crossbar is equal to the pass- 
band width of that detector). The filled circles are data presented by Newburn, Jr. 
and Gulkis (1971) . 

rv 1 

Fig. 9. Plot of k and G as functions of 05 and a. The ranges of the variables 
cover the ranges observed for the Jovian cloud layers . The values of k are 
identified on the right side of the graphs. For Jupiter k = 1, but values of k = 1.03 
and 0.97 are also given for comparison purposes. The values of G of each line 
is given within the graph . 

Fig. 10a - 10b. Plots of up versus wavelength and a versus wavelength, respectively, 
for the Jovian cloud layers. The open circles are for the STrZ and the filled 
circles are for the ER. 

Fig. 11 . Plot of k and G as functions of and a. The ranges of the variables 
cover the ranges observed for the Saturn cloud layers. The values of k are 
identified on the right side of the graph. The values of G of each line is given 
within the graph . 










' 0.5 1.0 1.5 2.0 


A _, ( jLim- 1 ) 
A (pm) 



0.5 1.0 1.5 2.0 

A"'(pnrr') 


i 

| 

i 

i 

i 

I 


FIGURE 4a-4b 










WAVELENGTH (pm) 



I 


WAVELENGTH (jum) 





K = 0.90 
K=0.85 

K=0.80 
K = 0.75 


FIGURE 11 



APPENDIX B 


THE SCATTERING MEAN FREE PATH 
IN THE 

URANIAN ATMOSPHERE 


* 


Michael J. Price 
Planetary Science Institute 
Tucson, Arizona 85704 


* Visiting Astronomer, Kitt Peak National Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under 
contract with the National Science Foundation 


- 2 - 



Page Intentionally Left Blank 



ABSTRACT 


New measurements of the equivalent widths of the 4-0 S(0) and S(l) H 2 
quadrupole lines in the Uranian spectrum have been obtained using high- 
dispersion (4.12A/mm) image-tube spectrography . The measured equivalent 
widths are 62-6 mA and 58-4 mA for the S(0) and S(l) lines, respectively. 

Curve-of-growth analysis in terms of a reflecting layer model yields an H 0 

+ 4-15 ^ 

column density of 780-200 km. amagat and a temperature of 78 jq °K. 

Interpretation using a semi-infinite, homogeneous, isotropically scattering 
model for line formation yields a scattering mean free path at A 6400A of 
550-150 km. amagat. The values derived for the H 2 column density and scatter- 
ing mean free path are insensitive to atmospheric temperature. In a clear, 
conservatively Rayleigh scattering, pure H 2 atmosphere, the corresponding 
scattering mean free path is 692 km. amagat. The results are consistent with 
the hypothesis that the regions of the Uranian atmosphere probed by visual 
radiation are essentially devoid of aerosol particles. 
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I. INTRODUCTION 


Some time ago Belton, McElroy and Price (1971) proposed that the region 
of the Uranian atmosphere probed by visible light is entirely devoid of clouds. 
Their proposition was based on an analysis of measurements by Giver and 
Spinrad (1966) of the 4-0 S(0) and S(l) E^ quadrupole lines. Interpretation in 
terms of both reflecting layer and simple scattering models for line formation 
indicated that the atmosphere can be represented by an essentially semi-infinite, 
clear, pure E^ model. Using a curve-of-growth analysis for the reflecting 
layer model, which included the effects of saturation, Belton et al derived an H 2 
column-density of 480 km. amagat and an atmospheric temperature of 118+ 40°K. 
In a detailed discussion of Raman scattering and pressure- induced dipole H 2 
opacity sources, Belton et al showed qualitatively that a semi-infinite, pure H 2 
model atmosphere is compatible with the general visible and near infrared 
spectral character of Uranus . To provide for the detailed band structure an 
additional contribution to atmospheric opacity by gaseous methane is required. 
Wallace (1972) has pointed out, however, that there are enough uncertainties 
in both the observational data and theoretical analyses to permit alternative views. 

Qualitative support for the clear atmosphere hypothesis comes from direct, 
high angular resolution (+0.2 arc sec.) imagery of the planet obtained by 
Danielson, Tomasko and Savage (1972) using a balloon-borne telescope. Except 
for the presence of strong, symmetrical, limb-darkening, the visual (X 3800A - 
X 5800A) appearance of Uranus is of a featureless disk. Although small 
quantitative discrepancies exist, the magnitude of the limb-darkening is consistent 
with predictions for a semi-infinite, homogeneous, non-conservative, Rayleigh 
scattering atmosphere. The planetary radius (25,900 kms) measured by 
Danielson et al is in good agreement with the optimum value (25,300 kms) derived 
by Belton et al from the ultraviolet geometrical albedo of Uranus on the basis 
of a semi-infinite, pure E^, Rayleigh/Raman scattering model atmosphere. 
Quantitative support is provided by new measurements by Lutz (1973) of the 
strengths of the 3-0 S(0) and S(l) H 2 quadrupole lines. Using curve-of-growth 
analysis for a reflecting- layer model, Lutz obtained several total line-of-sight 


column-densities v N(J), where i? is the geometrical "air-mass" factor and 
N(J) is the number of molecules with rotational quantum number J. Specifically, 
Lutz derived r? N(J=0) = 333+ 190 km. amagat and ?? N(J= 1) = 640+ 220 km. amagat. 
The corresponding atmospheric temperature was T = 111+30 °K. 

Nonetheless Prinn and Lewis (1973) have argued on thermodynamic grounds 
that a significant methane haze should be present high in the Uranian atmosphere, 
with a deep-dense ammonia cloud layer far below. They have also shown that 
the spectral character of the visible and near IR continuum could be explained 
in a straightforward way by a cloudy component in the atmosphere . On the 
basis of measurements of the infrared spectral albedo of Uranus, Binder and 
McCarthy (1972) suggested that a thin layer of CH^ crystals might indeed exist 
high in a predominantly clear atmosphere. Westphal (1972) has cited 
photoelectric-scanner evidence of limb-brightening in the strong methane band 
at 0.89 microns to support the hypothesis of an optically thin, high albedo, cloud 
layer in the Uranian atmosphere. Photographic confirmation of limb-brightening 
has been obtained by Sinton (1972). 

Limb-brightening is, however, predicted for the clear Hg-CH^ atmospheric 
model; it occurs in certain wavelength regions by virture of the vertical in- 
homogeneity which results from pressure-induced ppacity sources. Belton and 
Price (1973) have investigated in detail the variation with wavelength of the 
center-to-limb intensity profile for the case of a semi-infinite, non-conservative, 
inhomogeneous, Rayleigh scattering, - CH^ atmosphere. Measurement of 
the wavelengths at which the variation of intensity across the disk changes from 
limb-darkening to limb-brightening, and vice versa, provides a crucial test of 
the hypothesis that clouds are absent from the visible atmosphere of Uranus. 

All of the strong CH^ bands beyond 7600A should be limb-brightened, as should 
the central regions of the CH^ bands at 6190A, 6670A and 7300A. The observa- 
tions by Danielson et al , by Sinton, and by Westphal are not, however, sufficiently 
accurate and numerous, nor made in the critical spectral regions to be conclusive. 


To further explore the question regarding the clarity of the Uranian atmos- 
phere, the spectroscopic basis of the cloudless concept is re-examined in this 
paper. The 4-0 quadrupole band was selected for further study since the 
S(0) and S(l) lines lie in a spectral window between the A6190A and A6670A CH^ 
bands which is relatively uncrowded by Fraunhofer, or Telluric, line absorption. 
The 3-0 quadrupole band is less suitable since it is difficult to locate the 
continuum level in that spectral region (A ~ 8200A). New measurements of the 
equivalent widths of the 4-0 S(0) and S(l) quadrupole lines are reported. In 
this analysis, the scattering mean free path is indirectly determined. 


II. OBSERVATIONAL DATA 


Only Giver and Spinrad (1966) have published measurements of the 
strengths of the 4-0 S(0) and S(l) H 2 quadrupole lines in the Uranian spectrum. 
Using the Lick Observatory 120-inch aperture reflector/coude spectrograph 
system, they obtained two broadened spectra in the vicinity of A6400A at a 
reciprocal dispersion of 8A/mm. No information was given concerning spectrum 
width or spectral resolution. Individual uncertainties of +30 percent quoted for 
the derived strengths of both lines testify to the difficulty of the observations . 
Their measurements of the equivalent widths of both lines, as weighted by 
Belton et al (1971), are listed in Table I. 

Much of the problem in obtaining accurate equivalent widths for the H 2 
quadrupole lines results from their extreme narrowness. To improve the 
quality of the observational data, it was decided to obtain new broadened 
spectra at considerably higher dispersion using the maximum possible spectral 
resolution. Since high spectral resolution requires a narrow spectrograph 
slit, long exposure times are required. To maximize spectral resolution 
within a reasonable exposure time, it was decided to record the spectra using 
an image tube . 

During the night of 1972 May 22/23, an excellent quality high-resolution 
spectrogram of Uranus in the A6400* region was obtained with the Kitt Peak 
84- inch aperture reflector/coude spectrograph system. The spectrum was 
recorded using an ITT Type F4089 Image Tube located at the focus of the No. 6 
camera. Focussed with a permanent magnet assembly, this intensifier has a 
40 mm diameter photocathode on a plane parallel, ultraviolet transparent face- 
plate together with a phosphor output which is coupled by fiber optics to a 2 inch 
x 2 inch Kodak Ila-D plate. The maximum in the spectral intensity distribution 

o 

of the phosphor output lies in the region of A 5200A. During intensification, the 
original image on the photocathode does not suffer demagnification. Spatial 
resolution on the plate corresponds to 40 linepairs/mm (12.5 pm). Using the 
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first order spectrum, a reciprocal dispersion of 4.12A/mm ontheplate was 
achieved. During the course of the observations, the angular diameter of Uranus 
given by the 1972 American Ephemeris and Nautical Almanac was 3788 arc. 

During the exposure, the "seeing”, defined as the ha If -intensity diameter of a 
stellar image, varied in the range l"-2" arc. A slit width setting correspond- 
ing to an angular width of 0724 arc was chosen; the width projected on the plate was 

e 

38 pm. Spectral resolution on the plate was conservatively estimated to be 51 pm(0.21A) 
the sum of the spatial resolution of the image tube and the projected slit width. 
Broadening of the spectrum was carried out, using the No. 8 Decker, over an 
angular range of 476 arc; the corresponding width on the plate was 0.79 mm. 

A neon light source was used to provide a comparison spectrum. The spectro- 
gram of Uranus was exposed for 150 minutes. 

Calibration of the photocathode/phosphor of the image tube was carried out 
using the moon as a light source. To adequately record the lunar spectrum, an 
exposure of 30 seconds was required. Plate calibration was carried out with 
a standard spot sensitometer(13-step). By using a combination of intermediate- 
band color and neutral density filters, the spectral intensity distribution of the 

o 

spot sensitometer was matched as far as possible to the Uranian A6400A 
continuum produced by the phosphor of the image tube. Maxima in the respective 
spectral intensity distributions coincided within +1001 Exposure times were 
matched within a factor 3. Cut from the same original, the Uranian and lunar 
spectrograms were developed together with the spot calibration plate in D-76 
for 15 minutes at 68°F. In view of the wide disparity in exposure times for the 
Uranian and lunar spectrograms (30 secs-2 1/2 hours), prior experiments 
were carried out to determine whether the shapes and equivalent widths derived 
for the absorption lines might be affected by reciprocity failure in the phosphor / 
photographic plate combination. The experiments were based on recording the 
Uranian and lunar spectrum near A4000A. Using the Utrecht Photometric Atlas of 
the Solar Spectrum as a standard, no evidence of reciprocity failure was found 
throughout the relevant range in exposure times. 


/ 


- 9 - 


Density tracings were obtained using a Joyce-Loebl Automatic Recording 
Microdensitometer Mk.IIIC. For efficient handling of the data, the spectro- 
grams of both Uranus and the moon were digitized. Digital spacing was chosen 
to provide 2-3 data points per spectral resolution element. Following density- 
intensity conversion, Uranus/moon ratio spectra were obtained in the immediate 
vicinities of each of the H 2 quadrupole lines . Each spectrum was normalized 
to the continuum level associated with the center of the respective H 2 line. In 
each case, the normalization factor was determined using 30 adjacent continuum 
data points spread symmetrically on either side of the line . The runs error 
of their mean was taken as the uncertainty in the location of the continuum at 
the line center. As such, it may be used as a direct estimate of the uncertainty 
in the derived equivalent widths. Bear in mind, however, that only one plate 
was used in the determination of the equivalent widths. The observational data 
are presented in Figure 1 . Equivalent widths for both H 2 lines were determined 
using the trapezoidal rule, anchored at the digital data points. The results 
are listed in Table I. 
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HI. ANALYSIS 


Treatment of the problem of line formation in planetary atmosphere 
normally employs either of two distinct atmospheric models - the reflecting 
layer model or the scattering model. In the reflecting layer model, the 
atmosphere is assumed to consist of a clear layer of gas overlying a sharply 
defined, optically thick, cloud layer. The principal attraction of the model 
is that it utilizes a minimum number of parameters. Besides known line- 
strength data and geometrical factors, only the column-density of gas, the 
effective pressure, and the temperature are involved. Iterative application 
of curve -of -growth theory in the multi-line formation problem will permit 
the mean atmospheric temperature, effective pressure, and column-abundance 
of Hg to be derived. In the scattering model, the atmosphere is assumed to 
consist of a semi-infinite homogeneous, isotropically scattering, aerosol haze. 
Compared to the reflecting layer model, an additional parameter - the single 
scattering albedo in the continuum - must be defined; it is normally obtained 
by interpreting narrow-band photometric, or spectrophotometric, continuum 
measurements in terms of a semi- infinite, homogeneous, isotropically scatter- 
ing atmosphere. Also, the column-density of gas is replaced by the concept of 
the amount of gas in a scattering mean free path. Information concerning the 
mean volume density of aerosol particles above the effective level of photon 
penetration can be obtained. Both models will be used in the analysis. The 
Hg column-density obtained from the reflecting layer model will determine 
whether the effective depth of photon penetration is significantly reduced by 
“cloud particle extinction. The scattering model will provide an estimate of 
the average scattering mean free path in the regions of the atmosphere 
probed by visual radiation . 

Interpretation of the observational data is based on the curve-of-growth 
theory for homogeneous atmospheres, each model for line formation being 
individually considered. For both reflecting- layer and scattering cases, 
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single values of effective temperature and pressure are used to represent the 
Uranian atmosphere. The quadrupole lines' absorption profiles are 

determined not only by both Doppler-and pressure-broadening but also by the 
phenomena of collision-narrowing. Following Fink and Belton (1969) and Belton, 
Me Elroy and Price (1971), the Galatry (1961) profile will be used to describe 
the absorption profile. Revised values for the pressure-broadening coefficient 
at STP, o^, and the collision-narrowing parameter, z, were adopted from 
recent work by Belton (1973) on the 1-0 S(l) Hg quadrupole line. Computations 

by James (1969) indicate that the corresponding pressure-broadening parameter, 

-3 -1 -1 

a n , is 1.5 x 10 cm amagat . The value of is taken to vary with 

u -0 75 U 

pressure, p, and temperature, T, as pT ' . In the notation of Fink and 

Belton (1969), z is equal to 2ttcc^/P, where is the pressure-broadened 

Lorentz half-width (1/2 full line width), and P is the coefficient of dynamical 

friction. By definition p is equal to kT/mD, where k is the Boltzmann constant, 

m is the molecular mass, and D is the diffusion coefficient. On the basis 

of new information on the Hg diffusion coefficient obtained from Murray and 

Javan (1972), Belton derived a z value of 0.032. Using the revised values of 

cxq and z to compute a theoretical cur ve-of- growth, he satisfactorily fitted the 

laboratory data on the 1-0 S(l) line obtained by Rank, Fink and Wiggins (1966). 

Because of the lack of analogous information for the 4-0 S(0) and S(l) lines, the 

parameters derived for the 1-0 S(l) line will be used in this analysis. 

In the notation of Fink and Belton (1969), the relation between the E^ 
quadrupole line equivalent width, W, and the parameters of the radiative 
transfer model can be- written 


W 

tt JT) VF 


= F (u; P,^,y,z, £) 


(1) 


where is the Doppler width, and u is proportional to the characteristic optical 
length in the medium at the line center. In the isothermal reflecting layer model, 
the characteristic length is the scale height; in the scattering model it is the 
interaction mean free path. Specifically, u is equal to N S(J,T)/a D -/F, where 
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S(J,T) is the line "strength" containing the Boltzmann population factor for 
the lower level J of the transition at temperature T. In the reflecting layer 
model, the characteristic amount of gas, N, is the column-abundance; in the 
scattering model, it is the photon mean free path. The parameters p, Pq are 
the cosines of the angles of emergence and immergence, respectively; y is equal 
to a L / a D ; and w c is the single scattering albedo, not required in the reflecting 
layer model. 


Interpretation of the observational data in terms of the reflecting layer model 
closely followed the iterative curve-of-growth technique used by Fink and Belton 
(1969) to analyze H 2 quadrupole line data for Jupiter. Particular use was made 
of their "effective-pressure" concept. For the initial application of the curve-of- 
growth, an atmospheric temperature of 76°K was obtained on the assumption that 
the S (0) and S (1) lines were unsaturated; an effective pressure of 2 atm. was de- 
rived on the assumption that the reflecting layer was located at an H 2 Rayleigh 
scattering optical depth r (A.6400A) = 2. The relevant cross-section was derived 
from Victor, Browne and Dalgarno(1967) . Unit optical depth is equivalent to an 
H 2 column-density of 692 km amagat. Since the phase angle of Uranus was 
small, the values of p, Pq were taken equal at unity. Computations of S(J, T) 
were based on H 2 quadrupole absorption line data given by Dalgarno, Allison, 
and Browne (1969). With the atmospheric temperature held constant, iterative 
application of the curve-of-growth to the 4-0 S(l) line lead rapidly to an 
optimum effective pressure of 3.9 atm. , corresponding to a y value of 0.362. 
Iteration was based on the 4-0 S(l) line since the corresponding S (J, TJ/o^ -JiT 
factor is insensitive at the + 5 percent level to temperature variations throughout 
the range 75°K ^ T ^ 150°K; by comparison, the corresponding variation for the 
4-0 S(l) line is + 50 percent. The final curve-of-growth, together with the obser- 
vational data, is shown in Fig. 2. Taking into account saturation effects in both 


the S (0) and S (1) lines, the atmospheric temperature finally adopted was 

78 °K. The H 9 line-of-sight abundance derived from the 4-0 S(l) line was 

■10 +712 ^ 

then 2791 km. amagat. Quoted errors are based on the uncertainties in 


the measured equivalent widths only. 
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To derive the H 2 column-density, the line-of-sight abundance must be 
divided by the air mass factor, 77 . Because of seeing effects coupled with the 
small angular diameter of the planet, precise determination of 77 is not possible. 
Nevertheless, limitations may be placed on its value by considering two extreme 
cases. In both, the slit of the spectrograph crosses the diameter of the planet. 

In the first case, perfect seeing (i.e. no atmospheric angular spreading of near- 
point light sources) and symmetry with respect to position angle are assumed. 

The air mass factor is readily shown equal to 77 . The corresponding H 9 column 

+227 ^ 

density is 888 ^ km. amagat. In the other extreme, poor seeing (i.e. light 

from all regions of the planetary disk is completely "scrambled" before reaching 

the telescope) is assumed. The 77 value is 4. The resultant H 9 column-density 

+178 ^ 

is 698 ^ km. amagat. In fact, the seeing lies between the two extreme cases 

so that a mean 77 value is probably more nearly correct. The H 2 column-density 

finally adopted is 780+200 km. amagat. 

Application of the homogeneous, isotropically scattering model for line 
formation followed the general procedure used by Fink and Belton (1969) on the 
basis of theoretical developments by Belton (1968). Determination of the single 

/v 

scattering albedo, co c , was based on spectrophotometric measurements of the 
Uranian geometrical albedo, p, by Younkin (1970) corrected to the planetary 
radius obtained by Danielson et al (1972). For the case of a semi- infinite, 
homogeneous, isotropically scattering atmosphere, geometrical albedo is related 
to the single scattering albedo through the Chandrasekhar (1950) H-functions 
according to 

1 

P = ? “c|'b4 c <'V <*0 ( 2 ) 

The value derived for co was 0.975. Values selected for y and z were 0.362 and 

c . 

0.032, respectively. The atmospheric temperature was taken as 78°K. Because 
of insensitivity of the relevant S(J, T)/a^ -JT factor to uncertainty in the tempera- 
ture, determination of the scattering mean free path, N, was based on the measured 
equivalent width of the 4-0 S(l) line alone. 
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In analyzing the observational data, the two extremes of ’’seeing” were 
considered. For perfect seeing, the theoretical equivalent width, W, was 
obtained by averaging across the diameter of the planetary disk according to 



where 0 ( = cos p = cos * Mq) is both the angle of emmergence and immer- 
gence, W(0) is the corresponding equivalent width, and 1(0) is the corresponding 
specific continuum intensity. In the case of poor seeing, we have 



In each case, the relation between W and N was defined, trial values for the 
scattering mean free path of 1, 3, 10, 30, 100, 300, 1000, 3000, 10000 km. 
amagat being used in turn. Six -point Gaussian quadrature was used to perform 
the numerical integrations. The actual scattering mean free path (X6400A)_ 
was derived from the observational datum by graphical interpolation. For the 
case of perfect seeing, the scattering mean free path was 595 + 80 km. amagat. 
For the case of poor seeing, the corresponding value was 560+ 80 km. amagat. 
Quoted errors are based on uncertainty in the measured equivalent width. By 
carrying out similar computations for the extremes of temperature suggested 
by the reflecting layer model (i.e. 68°K and 93°K), a better estimate for the 
uncertainty in the scattering mean free path can be obtained. The value finally 
adopted was 550+ 150 km. amagat. Bearing in mind that in a pure conser- 
vatively Rayleigh scattering atmosphere the scattering mean free path 
(X6400A) is 692 km. amagat, neither model for line formation provides any 
evidence that aerosol particles inhibit the penetration of visual photons into 
the Uranian atmosphere. Results from the reflecting layer model suggest that 
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the deep-dense NH^ cloud layer predicted by Prinn and Lewis (1973) to occur at 
H 2 column-densities ~ 370 km. amagat must be much deeper, in fact too deep 
to be probed. Results from the scattering model support the hypothesis that 
the regions of the Uranian atmosphere probed by visual radiation are essentially 
devoid of cloud particles. If the aerosol haze predicted by Prinn and Lewis (1973) 
to occur at H 2 column-densities less than 119 km. amagat does in fact exist, it 
is physically and optically too thin to be detected by the Hg quadrupole lines' 
technique . 
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TABLE I 


h 2 QUADRUPOLE LINES 
MEASUREMENTS OF EQUIVALENT WIDTHS 


Line 

Wavelength 

(A) 

Equivalent Width (m A) 

Price 

Giver and Spinrad 

(4-0) S(0) 

6435.03 

62+6 

30+8 

(4-0) S(l) 

6367.80 

58+4 

41+ 15 

* 1 -9 
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